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Polarization transfers from an electron spin to a nuclear spin are essential for various physi¬ 
cal tasks, such as dynamic nuclear polarization in nuclear magnetic resonance and quantum state 
transformations on hybrid electron-nuclear spin systems. We present time-optimal schemes for 
electron-nuclear polarization transfers which improve on conventional approaches and will have 
wide applications. 
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I. INTRODUCTION 

As the gyromagnetic ratio of an electron is two to three 
orders of magnitude larger than the one of a nucleus, elec¬ 
tron spins are much easier polarized than nuclear spins. 
This offers a way to improve the polarization of nuclear 
spins by transferring polarization from electron spins to 
nuclear spins; much higher nuclear spin polarization can 
be achieved as compared to a direct polarization. This 
idea has been widely used in various physical settings, 
for example dynamic nuclear polarization (DNP) [IH5] 
employs this idea to dramatically improve the sensitivity 
of nuclear magnetic resonance (NMR) [II H] . It is also 
frequently used on various hybrid electron-nuclear spin 
systems, such as organic single crystals H], endohedral 
fullerenes [T(1f[T^ . phosphorous donors in silicon crystals 
[13j . and nitrogen-vacancy centers in diamond [HHn]. 
For example, in the case of nitrogen-vacancy centers in 
diamond, efficient polarization transfers are used to ini¬ 
tialize the quantum state of nuclear spins for quantum 
information processing. 

Efficient polarization transfers are practically achieved 
by properly engineered pulse sequences whose design is 
studied in the field of quantum control [I8l - l23] . In re¬ 
cent years, significant progress has been made in quan¬ 
tum control for both numerical and analytical 

[551 - HT] methods. Extensive knowledge has been gained 
on optimal pulse sequences for two- and three-level sys¬ 
tems [22HSS], two uncoupled spins [^[55], and two cou¬ 
pled spins [55H55] . Further advances have been made on 
how to optimally control multiple coupled spins [SSHSI]- 
These methods have been successfully applied in NMR 
[5^ 155] to designing broad-band and decoupling 

pulse sequences p7II102j . They have also been utilized 


in magnetic resonance imaging [25lll03Ill05] and electron 
paramagnetic resonance [TU5] . 

In this article, we consider time-optimal pulse se¬ 
quences for polarization transfers from an electron spin 
to a nuclear spin. Relaxation and decoherence are in 
practice inevitable and result in a loss of signal. But 
their effect can be mitigated by short pulse sequences 
which allow for highly sensitive experiments. We ana¬ 
lyze and explain how the form of time-optimal sequences 
depends on the direction of the polarization by studying 
time-optimal transfers for different directions. 

Recent analytical |107j and numerical [10811109] studies 
focused on low-field single-crystal experiments, where the 
nuclear Larmor frequency and pseudo-secular hyperfine 
interaction (see Sect. 3.5 of Ref. [T]) are comparable in 
magnitude. As in |110H113] . we focus here on the cases 
of secular hyperfine coupling (see Sect. 3.5 of Ref. m)- 
These assumptions are satisfied in liquid-state and high- 
held solid-state DNP. 

We analyze two particular cases of polarization trans¬ 
fers and determine the corresponding time-optimal se¬ 
quences. In Section]^ we consider the transfer from the 
state Sz of the electron spin to the state G of the nuclear 
spin. The second time-optimal transfer from Sz to 
is presented in Section |III| And most interestingly, the 
corresponding optimal transfer time is shorter by 78.5% 
when compared to the transfer from Sz to G, which high¬ 
lights that the transfer efficiency depends crucially on the 
target state of the nuclear spin. We discuss our results in 
Section IV and the possibility of a non-sinusoidal carrier 
wave form is entertained in Section We conclude in 
Section |V^ and certain details are relegated to Appen¬ 
dices El and m 
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II. TRANSFER FROM Sz TO h 


In this section, we study the polarization transfer from 
the initial state Sz to the final state Iz m- We assume 
a secular hyperfine coupling (see Sect. 3.5 of Ref. m)- In 
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the lab frame, the resulting Hamiltonian is given by 


H = UJsS^ + UJlh + ^TlASzh + + 2TTVa;{t)Ix, 

( 1 ) 

where ujs and toj denote the respective Larmor frequen¬ 
cies of the electron and the nuclear spin, A represents 
the strength of the secular hyperfine coupling, Ux{t) and 
Vx{t) are the amplitudes of the control fields. Here, Sj = 
(tjj (g)cro)/2 acts on the electron spin and Ik = {(jQ®ak)/2 
acts on the nuclear spin with j,k G {x,y,z}, where 
CTo := ( 0 1 ) denotes the identity matrix and the Pauli ma¬ 
trices are :=(%), (t^ := (° ~o), and := (J.?)- 
For typical NMR settings, only a single radio-frequency 
coil is used which can be assumed to be oriented along 
the X axis of the lab frame. Hence, only a single control 
Vx (t) appears for the nuclear spin in the lab frame Hamil¬ 
tonian of Eq. Q . We assume in this work that the carrier 
wave form for the nuclear spin has a sinusoidal shape, i.e. 
Vx{t) = v{t) cos[cLi}^t -I- with amplitude v{t) < 2uniax 
and phase (f>o. Here, cof is the carrier frequency of the 
radio-frequency irradiation and 2umax denotes the max¬ 
imal control amplitude (in the lab frame). This choice 
of Vx(t) is motivated by the properties (e.g., bandwidth 
limitations) of the usually available wave form genera¬ 
tors and amplifiers. More general carrier wave forms are 
discussed in Section IVl 

By switching to the rotating frame of ujsSz + 
corresponding to the carrier frequencies ujs and coy = 
oji — and applying the rotating wave approximation, 
we get an effective Hamiltonian 


Hzot — + Iz + 2ttASzIz 

+ + where (2) 

:= 2TrUx{t)Sx + 2-KUy{t)Sy, 

Hlit ■= ‘^'^Vx{t)Ix + 2-KVy{t)ly. 


One can obtain any desired offset term uj°^Iz in the drift 
term of i?rot in Eq. (p|) by suitably choosing the carrier 
frequency |115j . For simplicity, wj® is set to zero in 
the following. The microwave-frequency control pulses on 
the electron spin and the radio-frequency control pulses 
on the nuclear spin are given by and respec¬ 

tively. The control amplitudes Ux{t), Uy{t), Vx(t), and 
Vy{t) satisfy the bounds 


^ ^niax and + < 

^max; 

where Umax and Umax denote the maximal available am¬ 
plitudes of the control fields in the rotating frame for 
a given experiment. This is a result of the rotating 
wave approximation, which reduces the maximal con¬ 
trol amplitude of 2umax in lab frame to Umax in the ro¬ 
tating frame [8]. In the following, we will assume that 
Mmax ^ A ^ Umax and neglect the time needed to ap¬ 
ply operations that can be generated by the hyperfine 
coupling and the controls on the electron spin mg. 

Time-optimal transformations are essentially only lim¬ 
ited by the weak controls on the nuclear spin. The opti¬ 
mal strategy to achieve a desired transfer can be inferred 
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FIG. 1. Schematic depiction of the absorption profiles for 
(a) 2SzIz and (b) Iz- 



y 



FIG. 2. In the polarization transfer from 2SzIz to Iz, the 
P component of the nuclear spin doublet is rotated by an 
angle of tt around the y axis and the a component is left 
invariant (both visualized in the interaction frame). Note 
that the electron spin is in the state \p) on the left-hand side 
and in the state |a?) on the right-hand side. 


from the structure of cosets with respect to the fast op¬ 
erations [SniEgiUn]- Here, the fast operations are given 
by the hyperfine coupling and the strong controls on the 
electron spin. The transformation U which transfers Sz 
to Iz = USzU~^ will be suitably decomposed into a prod¬ 
uct U = U 2 U 1 . The unitary Ui transfers the initial state 
Sz to the intermediate state 2SzIz = UiSzUy^, and it 
can be generated using only fast operations. In addition, 
the unitary U 2 transfers the intermediate state 2SzIz to 
the final state Iz = U 2 { 2 SzIz)U 2 ^^ ^od one has to use 
the weak controls on the nuclear spin in order to gener¬ 
ate U 2 ■ Below, we will provide a time-optimal scheme to 
produce U 2 - This results also in a time-optimal scheme 
for U as any faster scheme for U would also imply a faster 
one for U 2 = UUy^. 

The polarization transfer from Sz to Iz can be decom¬ 
posed into the following steps mi- 


Sz 


-^Sx 


2SyIz — 


t 2SzIz Iz. (3) 


where we denote 


•S'" := (0 0) ® ^0 and 5 ^ := ([j ?) (g) cto, 

then irS^Iy = —nSzIy -|-7r/y/2. As shown in Eq. the 
polarization transfer from Sz to 2SzIz is accomplished 
using an INEPT-type transfer [51 lll8j : First, we apply a 
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hard 7r/2 pulse to the electron spin along the -\-y direc¬ 
tion (i.e. Sy). Then, we let the hyperfine coupling evolve 
for the duration of 1/(2^) units of time. Another hard 
7r/2 pulse on the electron spin along the +x direction 
completes the transfer to 2SzIz- All of these steps take 
negligible time, since they are either local operations on 
the electron spin or operations which can be generated 
by the coupling. In conclusion, we can completely focus 
on the last step in Eq. (|^ where we need to generate the 
propagator 

UP{0) = exp{-t9S^Iy) = exp[-i{-eSzIy + §4)] (4) 

for 6* = TT. The operator in the exponent of U^(9) in 
Eq. 0 is a single-transition operator mm- In partic¬ 
ular, the operator U^{-k) = exp{—iTrS^Iy) describes a 
transition-selective tt rotation around the y axis in the 
subspace spanned by the basis states |/3a) and |/3/3), 
where the subspace corresponds to the (3 component of 
the nuclear spin doublet ^ freqiwncy oji/{2Tr) -|- A/2 
[ng as shown in Figures and Here, |a) and |/3) 
are eigenstates of and 4, e.g., Sz\a) = \a)/2 and 
^zl/3) =-|/3)/2. 

In the following, we determine a time-optimal scheme 
to produce the unitary U^{'k). The set of all unitaries 
which transfer 244 to 4 are discussed in Appendix |A 1[ 
where we also show by extending the results in the cur¬ 
rent section that choosing a different element from this 
set of unitaries does not lead to a shorter transfer time. 

To determine the optimal transfer, we switch to the in¬ 
teracting frame of 2'kASzIz by applying the transforma¬ 
tion exp(i27rA44i)4ot exp(—i27rAS'z40- The Hamil¬ 
tonian of Eq. 0 changes to |120) 

H'lnt = + 2TrUx{t)[cos{TrAt)Sx — sin(7rAt)244] 

-I- 2TTUy{t)[cos{TrAt)Sy + sin(7rAt)244] 

-I- 2TTVx{t)[cos{nAt)Ix — sin(7rAt)244] 

-I- 2TTVy{t)[cos{nAt)Iy -I- sin(7rAf)244] (5) 

which can be also written as 

Hint = + 2TTUx{t)[cos{TrAt)Sx — sin(7rAt)244] 

-I- 2TTUy{t)[cos{TrAt)Sy + sin(7rAt)244] 

-I- 2TT[vx(t) cos(7rAt)4 + Vy{t) sin(7rAf)244] 

— 2Tr[vx{t) sin(7rAf) — Vy(t) cos(7rAt)]5'“4 

-I- 2Tr[vx{t) sin(7rAf) -f Vy(t) cos^irAt)]S^ly, (6) 

where S'^Iy = Szly + 4/2 and S^Iy = —Szly + Iy/2. 

We aim at generating the operator in minimum 

time, which corresponds to maximizing the coefficient 
2TY[vx{t) sin(7rAt) -f Vy{t) cos(7rAt)] in front of S^Iy. The 
Cauchy-Schwarz inequality implies 

[vxit) sin(7rAt) -|- Vy{t) cos(7rAt)]^ < 

[vlit) + vl{t)][sm'^iTrAt) + cos{7rAt)^] < 


where the second inequality is a consequence of the con¬ 
straints on the amplitude of the control fields. The max¬ 
imal value of 27r[ti2,(t) sin(7rAt) -f Vy{t) cos{ttA t)] is de¬ 
noted by 27runiax and it can be achieved by choosing 
the controls Ux(t) = Uy{t) = 0, Vx(t) = Umaxsin(7rAt), 
Vy{t) = Umax cos (tt At). To Understand that this choice 
generates the desired operator, one can substitute the 
controls in the Hamiltonian with the chosen values and 
obtains 

i/int = -b 27rUmax sin(7rAt) cos(7rAt)4 

-b 27rumax cos(7rAt) sin(7rAt)2S'z4 
+ STTUmax COs{2t:A t)S°‘Iy + 27rUmaxS'^4- 

Since A ^ Umaxj average Hamiltonian theory implies that 
the first three terms average out to zero; and one is left 
with the desired Hamiltonian 27rUmax<S'^4- 

The minimum time to generate U^{tt) is then fixed by 
the relation 27rUmax4nin = tt, and one obtains 


The presented time-optimal control corresponds to a 
radio-frequency irradiation at frequency uij/{2tt) + A/2 
with duration Tminj which results in a transition-selective 
inversion of the /3 line of the nuclear spin doublet. This 
belongs to the class of controls presented in Ref. m 
and is also closely related to selective population inver¬ 
sion (SPI) experiments |121H124] . 



5 10 15 20 5 10 15 20 

Time (jjs) Time (jjs) 


FIG. 3. Numerically optimized pulses for the polarization 
transfer from Sz to 4. The coupling strength is 10 MHz and 
the bounds on the micro-wave and radio-frequency amplitudes 
are Umax = 1 MHz and Umax = 20 kHz, respectively. The 
maximal transfer efficiency is reached after 25fis. The insets 
show magnified parts of the controls Vx{t) and Vy{t) in order 
to illustrate their form. 
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We can also compute the maximal transfer efficiency 
Vmax{T) for a given time T. The operator U^{ 9 ) transfers 
the state 23^1 z to the state 

Uy^{e){2Sziz)[u^io)]^ = 
cos^{^)2SzIz + cos(|) sm{^)2SzIx 
-cos(f)sin(|)4+sin^(|)/^. (7) 

For 0 = 2'KVyna.xT, we get the maximal transfer efficiency 

VmaxiT) = sin^d) = sin^(7r'Umax'r) (8) 

for the transfer to Iz- Note that rjmaxiTmin) = 1- 

We compare our analytic results with numerical op¬ 
timizations for achieving the transfer from Sz to Iz as 
shown in Fig. § cf. nuHnsi- For these optimizations, 
the hyperfine coupling constant is chosen as ^ = 10 MHz 
and the maximal allowed radiation amplitudes are set to 
Umax = 1 MHz and Umax = 20 kHz |125j . In Fig. the 
transfer is completed after 25/is = l/(2umax) units of 
time which agrees with the analytically computed time. 
Moreover, the form of the numerically optimized controls 
compares nicely with the analytic results: the values of 
Ux(t) and Uy(t) are most of the time small (except for the 
beginning), and Vx(t) and Vy(t) have a sinusoidal form 
with the maximal allowed amplitude. 


III. TRANSFER FROM Sz TO OR 4 


We analyze now how to time-optimally transfer po¬ 
larization from the state Sz to 4 (and similarly for the 
transfer to ly). The considered transfer consists of the 
following steps: 






t 2SyIz 


A 2SzIz 


TL qa T _‘JL q& T 

2 >-> Jy 2 


t Ix 


where ^S°‘Iy — ^S^Iy = irSzIy As in Sec. we 
can focus on generating the final propagator U 2 = 
e:iqy{—iTrSzIy) in the product U = U 2 U 1 . The transfer 
Ix = USzU~^ is decomposed into a fast transfer to the 
intermediate state 2SzIz = UiSzlI^^ and a slow transfer 
to final state 4 = U2{2SzIz)U2^ ■ Building on the re¬ 
sults in this section, we prove in Appendix |A 2| that one 
cannot reduce the transfer time by substituting U 2 with 
a different unitary V satisfying 4 = V{ 2 SzIz)y~^ ■ 
Previously, the propagator exp(—zTrS'j/y) in the final 
step has been achieved [nni by applying a transition- 
selective radio-frequency —7r/2 pulse along the y direc¬ 
tion at the /? transition with frequency ujj/{2tt) -I-A/2 as 
well as a transition-selective radio-frequency 7r/2 pulse 
along the y direction at the a transition with frequency 
wi/{2tt) — A/2, see Fig. In the rotating frame of 
Eq. this irradiation scheme on the nuclear spin cor¬ 
responds to a radio-frequency Hamiltonian of the form 


= - Stt- 


■ cos(7rAt)4 — 27r AIAii sin(7rAt)4 


/3 ^ a. 




FIG. 4. In the polarization transfer from 244 to Ix the /3 
component of the nuclear spin doublet is rotated by 
around the y axis and the a component is rotated by 7r/2 
around the y axis (both visualized in the interaction frame). 


-\- 27rcos(—7rAt)4 + 27rsin(—7rAt)4 

= - 27rUniax sin(7rAt)4. 

Note that in this scheme the a and /? transitions can 
only be irradiated with a radio-frequency amplitude of 
'Cmax/2 in order not to exceed the maximal available 
radio-frequency amplitude Umax for the overall irradia¬ 
tion at the nuclear spin. Hence, the duration for the 
simultaneous ±7r/2 pulses along the y direction at the 
a and P transitions is equal to l/(2umax)- This conven¬ 
tional transfer is optimal if one considers only pulses at 
the frequencies a;//(27r)±A/2 of the nuclear-spin doublet 
(refer to |linj and the discussion in Section [rv| . 

Here, we show that shorter pulses are possible if one 
considers more general transfer schemes. Without ex¬ 
ceeding Umax, shorter pulses can be obtained by irradiat¬ 
ing at the frequencies uji/{2tt) ± A/2 with higher inten¬ 
sity since the resulting higher amplitude can be then de¬ 
creased by irradiating at additional well selected frequen¬ 
cies. Our approach is quite effective although it might 
seem counterintuitive at first. 

In the interaction frame of 2'kASzIz-, the Hamiltonian 
is again given by Eq. ([^. In order to generate the op¬ 
erator exp^—iirSzIy) in minimum time, we maximize the 
coefficient —2TrVx{t) sm{7rAt) of 2SzIy Note that 

—2TrVx{t) sm{TrAt) < 

\-27rVxit) sin(7rAt)| < 27rumax| sin(7rAt)|, 

where the second equality is implied by the constraint 
|ua;(t)| < Ux{t) + Vy{t) < Umax On the control ampli¬ 
tudes. Therefore, the maximal value 27rumax| sin(7rAt)| 
for —27rua;(t) sin(7rAt) can be attained by choosing the 
controls Ux(t) = Uy(t) = Vy(t) = 0 and Vxit) = 
—sgn[sin(7rAt)]umax- This means that Vx{t) is a square 
wave such that Vx(t) = Umax when sin(7rAt) < 0 and 
w(f) = ~^^max when sin(7rAt) > 0. In this case, the 
radio-frequency Hamiltonian in the rotating frame is 
given by 


= -27rUmax sgn[sin(7rAt)] 4 


2 


2 


(9) 
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FIG. 5. Numerically optimized pulses for the polarization 
transfer from to l^. The coupling strength is 10 MHz and 
the bounds on the micro-wave and radio-frequency amplitudes 
are given by Umax = 1 MHz and Umax = 20 kHz, respectively. 
The maximal transfer efficiency is reached after 20/iS. The 
inset shows a magnified part of the control Vx{t) in order to 
illustrate its form. 


We obtain the minimum time 

Tmin ~ / (St^max) 

for generating expl—inSzIy) in the interaction frame. 
The duration of the transfer is reduced to 78.5% of the 
length of the conventional pulse sequence. By transform¬ 
ing the operator back to the rotating frame, we obtain the 
operator exp(—i(j)SzIz) s^p{—iT^SzIy) e^p{i(l>SzIz) where 
(j) = 27r2irniin denotes the phase accumulated during the 
time Tmin = 7r/(8uniax)- The effect of this superfluous 
phase (p can be reversed using the hyperhne coupling 
2ttASzIz which takes only a negligible time period as 
the coupling strength A is much larger than the control 
strength Umax of the nuclear spin. Thus, the minimum 
time in the rotating frame is also given by 7r/(8uinax)- 
Similarly as in Eq. Q, we compute the maximal transfer 
efficiency 

Vmax{T) = sin(4UmaxT) (10) 

that can be reached for the polarization transfer from Sz 
to Ix in a specified time T. 

Transferring the state from Sz to ly is similar. We set 
Ux{t) = Uy{t) = Vx{t) — 0 and maximize the coefficient 
—2'KVy(t) sin(7rHt) of —2SzIx in Eq. ([^ by setting Vy{t) = 
—sgn[sin(7rHt)]uinax- The minimum time for this case is 
also given by tt/( 8uinax)- 

A numerically optimized pulse sequence for transfer¬ 
ring polarization from Sz to Ix is shown in Fig. cf. 


[mHis]. The maximal transfer efficiency is reached af¬ 
ter 20/is which is consistent with the analytical result of 
7r/(8Uinax) « 19.635^5. 

IV. DISCUSSION 

We see that the minimum time for transferring Sz to 
Ix or ly is shorter by a factor of 7r/4 when compared to 
the minimum time for transferring Sz to Iz- This factor 
can be explained by a closer examination of the pulse 
sequences. The radio-frequency sequence for the trans¬ 
fer from Sz to Iz shows a sine-cosine wave modulation of 
maximal amplitude for the Vx- and Uj,-components (see 
Fig. §. However, the radio-frequency sequence for the 
transfer from Sz to Ix consists of a square wave of max¬ 
imal amplitude for the U 2 :-component of the control (see 
Fig.§. The higher effective amplitude at the two fre¬ 
quencies uji/{2tt) ± A/2 and the shorter transfer time 
can be explained by decomposing the square wave into a 
sum of sine waves: 

/square(t) = sgn[sin(7rAt)] = - - sin(nAt). 

TT n 

n odd, n>l 

This is illustrated in Fig. where the first sine wave 
function has an amplitude which is larger by a factor 
of 4/7r when compared to the amplitude of the square 
wave. Therefore, the square wave contains implicitly a 
sine wave with an higher effective amplitude. This im¬ 
plies that the duration of the simultaneous ±7r/2 rota¬ 
tions of the a and /3 components of the nuclear spin dou¬ 
blet (see Fig. is shorter by a factor of 7r/4. 



Time 

FIG. 6. Decomposition of a square wave into sine waves: a 
large number of harmonics sum to an approximate square 
wave. The first harmonic has an amplitude of d/vr while the 
square wave has an amplitude of 1. 
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FIG. 7. The maximal transfer efficiency (i.e. fidelity) is shown 
for the transfers from Sz to Iz and 1^ in (a). The conpling 
strength is 10 MHz, and the control strengths are Wmax = 
20 kHz and Vmax = 1 MHz. In (b), data points for the transfer 
from Sz to Iz are shown also for the bandwidth-limited cases 
with bounds of 5.0 MHz, 4.8 MHz, and 4.6 MHz. Similarly, 
data points for the transfer from Sz to Ix are shown in (c). 


The square-modulated transfer sequence is optimal but 
needs infinite bandwidth. We also studied numerically 
how the maximal transfer efficiency varies as a function 
of time and bandwidth limitations. The results are shown 
in Fig.[7[ cf. |111H113] . In the case of infinite bandwidth, 
the results are consistent with the analytical results. The 
transfer functions sin^(7r'(;niaxi) and sin(4umaxi) for the 
respective transfers from Sz to Iz and Ix have been ob¬ 
tained in Eqs. Q and ( [To| . 

We compare our results to the time-optimal synthe¬ 
sis of unitary transformations in m- Motivated by 
energy considerations, only irradiations at the two res¬ 
onance frequencies a;//(27r) ± yl/2 of the control system 
were considered in |110j . This did not allow for the faster 
scheme obtained in Section IIIII which has been also ob¬ 
served numerically in |111H113| . The numerical results 
in Fig. also show that the faster scheme of Section |III| 
is—in a strict sense—only applicable in the case of infi¬ 



FIG. 8. (Color online) Normalized amplitudes of the short- 
time Fourier transform for different controls: (a) u-controls of 
Fig. m (b) w-controls of Fig. (c) u-controls of Fig. and (d) 
u-controls of Fig. 1^ In (c), (essentially) only the characteristic 
frequency of A/2 = 5 MHz is present. This is in contrast to 
(d) where also multiples of the frequencies ±A/2 appear. 


nite bandwidth. It provides in general superior results, 
but its benefit depends on the available bandwidth. 

The irradiations at the different frequencies can be 
clearly observed in Fig. where the normalized ampli¬ 
tudes of the short-time Fourier transform |126] are plot¬ 
ted for the relevant cases (using the method and imple¬ 
mentation of |127j ). The important difference between 
the controls of Fig. for the polarization transfer from 
Sz to Iz and the controls of Fig. for the transfer from 
Sz to Ix manifest itself in the short-time Fourier trans¬ 
forms of the V parts of the controls which are visible in 
subfigures (c) and (d) of Fig. In subfigure (c), one 
notices the characteristic frequency of A/2 = 5 MHz. In 
contrast to sub figure (c), many more frequencies appear 
in (d) at multiples of the frequencies ±A/2. This agrees 
with the square-wave form of the control Vx in Fig. 

In general context, our results can also interpreted as a 
connection between energy and time optimizations. The 
energy optimization leads to a sinusoidal solution while 
the time optimization leads to a square wave (see Fig.[^. 
This phenomenon has been also observed for the simul¬ 
taneous inversion of two uncoupled spins [571 [5H] where 
the minimum energy solution was related to the first har¬ 
monic in the Fourier expansion of the time-optimal solu¬ 
tion. The same reasoning applies here to the solutions of 
Section nn 
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We discuss now the possibility and opportunities of 
non-sinusoidal carrier wave forms for the electron and 
nuclear spin. Here, we focus on the nuclear spin as the 
maximal amplitude fmax limits the minimum polariza¬ 
tion transfer times in the rotating frame. But similar 
arguments might also be used for the micro-wave carrier 
wave form in applications where the minimum duration 
of an experiment is limited by Umax- As explained in 
Section [TTj we have considered so far a sinusoidal carrier 
wave form Vx{t) = v{t) cos[a;j^t-|-(^(t)], which is motivated 
by the limited bandwidth of typical radio-frequency wave 
form generators and amplifiers. Equation (|^ states the 
time-optimal radio-frequency Hamiltonian for transfer¬ 
ring polarization from the state Sz to Ix in the rotat¬ 
ing frame. In the lab frame, the corresponding radio¬ 
frequency Hamiltonian is given by 

= dTTUniax cos{ujft — 7r/2) sgn[sin(7rAt)] lx- (11) 

However, it is conceivable (e.g., for applications at low 
magnetic fields or for nuclei with small gyromagnetic ra¬ 
tios) that the resonance frequency and the correspond¬ 
ing carrier frequency of the controls are sufficiently 
small such that non-sinusoidal wave forms (containing 
higher harmonics of the carrier frequency ujf) can be cre¬ 
ated and amplified. One can therefore envision a radio¬ 
frequency Hamiltonian for the ideal case of infinite band¬ 
width as given by 

= dTT'Cmax sgn[cos(w}^t-7r/2)] sgn[sin(7rAt)] 4 . (12) 


Assuming the same maximal radio-frequency 2vjna,x (in 
the lab frame) and switching from the Hamiltonian 
in Eq. 0 to the Hamiltonian 77"'^ in Eq. ( |I^ , the radio¬ 
frequency amplitude of the carrier frequency is implicitly 
increased in the rotating frame by another factor of 4/7r 
(similarly as in discussed in Section IV). Consequently, 
the polarization transfer from the state Sz to 4 would be 
achievable using only (7r/4)^ « 61.7% of the conventional 
transfer time. 


VI. CONCLUSION 
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Appendix A: Decomposition of unitaries 
1. Unitaries which transfer 2SzIz to 4 


All unitaries in SU(4) can be be decomposed as K 1 AK 2 
with a slow evolution A := exp[—i(aS°‘Iy -+- /3S'^4)] and 
fast unitaries Ki and IC 2 which can be generated by con¬ 
trols on the electron spin and the secular hyperfine cou¬ 
pling (cf. [Ml ESI 111011128] !: recall that 5“ = (J 0 ) G erg 
and = (0 1 )G)cro. The unitaries that transfer 2SzIz to 
4 can be determined as solutions to the matrix equation 
4 = KiAK2(2SzIz)R^2^^■ The fast unitaries Ki and 
-K 2 can be parameterized using canonical coordinates of 
the second kind (see Section 2.8 of |129j . Section 2.10 of 
[130j . or Ghapter HI, Section 4.3 of |131j 'l. i.e. 


ATi := e 


icLjSx^ ^ ici^Sz 

3 ‘iQy4‘^S X IZ ^ id^'ZSy I z z 1 Z c> Z 


(Ala) 


AT, :=e 


— ibiSz „—ib2Sy —ib32Sz:Iz „—ibi2SyIz 


X g-ib^Sz ^-ibfi2SzIz ^-ibjlz 


(Alb) 


The surjectivity of the representations in Eq. (A1) is veri¬ 
fied in Appendix[^ As the unitary Ki commutes with 4 
and parts of K 2 commute with 2SzIz, tbe matrix equa¬ 
tion simplifies to 


With the help of the computer algebra system Maple 
[132j . one can verify that either a = 2 ttzi and /3 = tt -|- 
2t:z2 or a = 7r-|-27rzi and j3 = 27:Z 2 with zi,Z 2 G Z holds. 
In Sec. |ll]of the main text, we focused on the first case 
assuming that a = 0 and /3 = tt (i.e. zi = Z 2 = 0), all 
other cases are similar. 


We have presented time-optimal polarization transfers 
from an electron spin to a nuclear spin for the case of 
secular hyperfine couplings. In particular, we have ana¬ 
lyzed the transfers from the electron-spin state Sz to the 
nuclear-spin states 4 and 4- For the transfer to 4, we 
could improve on the duration of on-resonance sinusoidal 
solutions by applying a control which has the form of a 
square wave. Our results also highlight differences be¬ 
tween optimizations for minimum energy and minimum 
time. We have also discussed how these differences are 
related to bandwidth limitations. 


2. Unitaries which transfer 2SzIz to 4 

Similarly as in Appendix |A 1[ the unitaries can be 
decomposed into a product K 1 AK 2 of fast unitaries 
ATi, K 2 and a slow evolution A = exp[—i{aS°‘Iy -\- 
l3S^Iy)]. In particular, all unitaries which transfer 
2SzIz to 4 have to satisfy the matrix equation 4 = 
KiAK 2 { 2 SzIz)K 2 ^^^ 1 - observing trivial commuta¬ 

tors, the matrix equation simplifies to 

J _ g-iai2Sz:Iz ^-ia.32SyIz ^-ia22SzIz g-iailz 













^ ^-ibiS:^ ^-ib2Sy ^-ib32SxIz ^-ib42SyIz f2S I ) 

X ^'^b^2SxIz g^^2‘S'y j4^ 

X 2S^ Iz g*fi32Sy ^icL4^2SxIz 

With the help of the computer algebra system Maple 
[132) . one can infer that /? = a — tt + 27r2: holds for z ^TL. 
In Sec. [m] of the main text, we consider the case of A = 
exp[— 7 rzS' 2 /y] which corresponds to a = 7 r/ 2 , (3 = — 7 r/ 2 , 
and z = 0. This choice is actually optimal: It follows 
from Eq. (rol) in the main text that 


_ g iO'A‘ZS^I z ^ z 

^ g iOib ^z g Iz g^*^5 ^z g id’bSz g io,‘2 ‘2Sz Iz g^^5 Sz 

_ ^—iaA[cos{ab)2Sz;Iz+sin{ab)2SyIz\ 

X g-“3[cos(a5)2Sj,/^-sin(a5)2Sj,/j]g-ia22Sj/^ 

_ Q-ia'l 2SzIz g“®“3 2.SyIz g-»“2 25^ 


where the last step follows from the Euler-angle decom¬ 
position. Similar arguments for the conjugations with 


g lasSy g zaySi demonstrate Eq. (Bl). Any ele- 


ri ment in the connected Lie group that is infinitesimally 

= / —27r[u2,(t) sin( 7 rAt) — cos( 7 rAt)] dt, generated by the elements —iS^, —iSy, —iS^, —i2S, 

/)OC^rr«na.iOC*r/^«*>nl-»/~>r» <-iV» t /Ai VvT T r\ r\-r 


13= 27r[ua;(t) sin(7rAt) -|- Vy{t) cos(7rAt)] dt 

Jo 

holds for any given time T. Consequently, (3 — a = 
iTTVx{t)sm{TTAt)dt. One applies the condition /3 = 

a — n + 2ttz and obtains fj' 47rVx(t) sin(TrAt)dt = — tt -|- 
27r2:. This implies that 

I f 4:TTVx{t) sin^nAt)dt\ > | — tt -|- 27rz| > tt. 

Jo 

On the other hand, one has \ 4TTVx(t) sui(TTAt)dt\ < 
sin( 7 rAt)|dt = Sv^naxT (see Sec. 


Ill I. In order 


fo^ 47ru, 

to satisfy the condition /3 = a — tt + 27rz, the time T has 
to fulfill the inequality Sv^axT > tt. One gets a lower 
bound Tmiii > 7r/(8ui„ax) on the minimum time Tmin. In 


summary, the scheme presented in Sec. |III| of the main 
text is optimal as it saturates the lower bound. 


Appendix B: Verification of the surjectivity of the 
representations in Eq. (All 


In order to verify the surjectivity of Ki in Eq. (Ala) 
it is sufficient to verify the surjectivity of the product 

All = ^ 1 ( 07 , 06 , 05 , 04 , 03 , 02 ) := 

g ZCI 7 Sx g g id^Sz g i(i^2SxIz g icL^2SyIz g ici2‘2iSz^z 


—i2SyIz, and —i2SzIz can be achieved by a finite prod¬ 
uct of elements having the form of Ki\ this is a con¬ 
sequence of Lemma 6.2 in |133j . We apply Eq. (Bl) 


and the Euler-angle decomposition multiple times and 
obtain Ari(o 7 , 05 , 05 , 04 , 03 , 02 )Ali(d 7 , 05 , 05 , 04 , 03 , 02 ) = 
Ari(c 7 , Ce, C 5 , C 4 , C 3 , C 2 ) for certain values of C 7 , Ce, C 5 , C 4 , 
C 3 , and C 2 . In summary, we have verified the surjectivity 
of the representations Ki and Ki. 

Similar as for Eq. (Bl), one can verify that 


g ici^Sz^ icL^2SyIz 


holds for some 04 and O 3 . Consequently, the surjectivity 
of Ki implies the surjectivity of K 2 . 

An alternative second argument for the surjectiv¬ 
ity of Eq. (Ala]) applies the decomposition K'^A'Kl^ 


for the set K = exp(£) of all fast operations where 
AT' = exp(f') and A' = exp(a'). This decompo¬ 
sition is a consequence of the Cartan decomposition 
i = t' 0 p' where the corresponding linear subspaces 
are given by f' := span{— iSx, —iSy, —iSz, —Hz}, P^ := 
span}—i2SxIz, —i2SyIz, —i2SzIz}, and the abelian sub¬ 
algebra a' := span{—i25'2/z} C p' |128) . The decomposi¬ 
tion K[A'K 2 implies that the decomposition 

JJ> _ jj^inSzIz _ ^-idiSz! ^-id2Sy ^-idbSz Q-idi2SzIz 
X g id^Sz^ ^df^Sx^ '^djSy^ id^dz 


which consists of the first six elements of ATi as the sev¬ 
enth element commutes with all the other ones. First we 
show that there exists 04 , Cg, and such that 

g Sx g io>G Sy g idbSz g idA2Sx dz g icLb2Sydz g io>22Sz dz 

g ^^,^23 X d z g iQi,^2S y d z g iQ>22S z d z g idj Sx g Sy g idb S z 

(Bl) 

which can be written as 

g i’dA2Sxdz^ ^d,A‘2Sydz^ i(l22Szdz _ g idjSx^ ^dbSy ^ idbSz 

X (g“*“4 2Sx/2 ^-id32Sydz ^-id22Szdz'^^idbSz ^idsSy ^idjSx 

The effect of the conjugation with exp(— 1055 ^) is 

^—idbSz ^—idi2Sxdz ^—id32Sydz ^—id22Szdz ^idbSz 


is a surjective parameterization of the set of all fast op¬ 
erations. Therefore, the surjectivity is also verified for 

JJ _ ^-idiSx ^-id2Sy ^-idbSz ^-idi2Szdz 

X g ^db Sz g idg Sx g idjSy g id,^ dz g ittSz dz 

_ ^—idiSx g—id2Sy ^—i{d3+db)Sz ^—i{d4+ir/2)2Szdz 
^ ^ide2Sydz ^-idj2Sxdz ^-idgdz 

_ g id±Sx g id2Sy g 'i'd^Sz g ^dA2Sxdz 

X ^~^d 32 Sydz ^—id 32 Szdz ^—^dsdz 


where the last equality follows from the Euler-angle de¬ 
composition. This completes the second argument for 
the surjectivity of Eq. (Ala). 
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